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Abstract 

Improving the best known examples, two planar straight-line graphs which cause the non-termination 
of Ruppert's algorithm for a minimum angle threshold as low as a g; 29.06° are given. 

Introduction 

Given a planar straight-line graph (PSLG), Ruppert's algorithm [5] produces a conforming Delaunay tri- 
angulation satisfying a minimum angle bound. The standard analysis jH [10] demonstrates that when the 
input contains no angles smaller than 60° Ruppert's algorithm produces a size-optimal mesh (up to a con- 
stant factor) for any minimum angle bound a ;$ 20.7°. A more detailed analysis can slightly improve this 
restriction to a J$ 22.2° for non- acute input [5] and an additional (very mild) assumption further improves 
the guarantee to 26.5° [3]. 

In practice, Ruppert's algorithm succeeds for substantially larger minimum angle bounds than those 
guaranteed by the theory. Ruppert observed that the minimum angle reaches 30° during typical runs of the 
algorithm 8 . In a number of experiments Shewchuk found the algorithm to terminate for a Js 33.8° [5]. 
Certain modifications of the vertex insertion procedure further improve the algorithm suggesting that this 
constraint can be improved to possibly 40° or more [T]. 

Pav demonstrated an example of non-termination of Ruppert's algorithm on a simple non-acute PSLG 
for any a > 30° [4]. This example combined with the fact that the analysis of Ruppert's algorithm on point 
sets breaks down at 30° led to a natural conjecture that Ruppert's algorithm terminates and produces a 
well-graded mesh for all a < 30°. While some examples corroborated this idea [6], a recent example shows 
non-termination for a as low as « 29.51° [7]. 

We give two examples that improve upon those given in [7] . The first example produces non-termination 
of Ruppert's algorithm for a gl 29.10° and has a minimum input angle of about 87°. The second example 
gives a slight improvement, a <J 29.06° but requires a 60° input angle. 

First, some notation is defined. The line-segment between endpoints V! and v 2 is denoted viv 2 and 
the triangle with vertices Vi, v 2 , and V3 is Aviv 2 V3. Let Zviv v 2 be the angle at vertex vo between 
line-segments vTyvT and vqv 2 . 



Example 1 

This example PSLG consists of five input vertices and three adjacent input segments. By carefully construct- 
ing the input (as described below) Ruppert's algorithm inserts a circumcenter, followed by three consecutive 
circumcenters that yield to midpoints of the segments. The result is a configuration that is similar to the 
input but exactly half the size and thus Ruppert's algorithm can repeat this cycle indefinitely. Depicted in 
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Figure 1: Example 1. The black vertices and segments comprise the input. Blue vertices are circumcenters 
of skinny triangles: the dark blue vertex is inserted into the triangulation while the light blue vertices are 
rejected for encroaching segments. The red vertices are midpoints of encroached segments which are inserted 
by Ruppert's algorithm, (left) The sequence of vertices considered for insertion by Ruppert's algorithm with 
the relevant skinny triangles and encroached diametral circles, (right) Labels for each of the vertices. 

Figure [1] the sequence involves four vertex insertions: (1) the circumcenter ci of skinny triangle AV0V1V4 
is inserted, (2) the circumcenter C2 of skinny triangle AV0V2C1 encroaches segment V0V2 which causes the 
insertion of midpoint mi, (3) the circumcenter C3 of skinny triangle AV0V31111 encroaches segment vovi" 
which causes the insertion of midpoint ni2, and finally (4) the circumcenter of skinny triangle AvoVim2 
encroaches segment vovT which causes the insertion of midpoint ni3 . 

Construction The following construction describes how to place the input vertices (vo, . . . , V4) to achieve 
non-termination for the smallest possible angle threshold. We begin by fixing 71 which serves as the smallest 
angle for a number of skinny triangles in the example. Then the placement of the vertices is determined by 
the following steps; recall Figure [TJ 

1. Begin with a segment between two vertices Vo and vi. 

2. Let 013 be the midpoint of vovi. 

3. Place v 4 such that Zv4V 1 v = 90° and ZV1V4V0 = 71. 

4. Let Ci denote the circumcenter of AV0V1V4. 

5. Select V2 such that the circumcenter of AV0V2C1, denoted C2, lies on the diametral ball of voV2- 

6. Let mi denote the midpoint of vovj- 

7. Select 1x12 so that the circumcenter of AvoViiri2, denoted C4, lies on the diametral ball of vovi and 
the opposite side of vovT as ci. 

8. Define v 3 so that m 2 is the midpoint of V0V3. 

9. Let c 3 denote the circumcenter of AV0V3III1. 
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When a > 71 the construction ensures that triangles AvoVi V4 and (following the insertion of ci) AV0V2C1 
are split by Ruppert's algorithm. Moreover, if ni2 is inserted, then AvoVini2 is skinny. Also, by construction, 
if/ when C2 and C4 are inserted, they encroach V0V2 and vovT, respectively, resulting in the insertion of mi 
and m 3 . 

For non-termination to occur, AV0V3III1 must also be skinny and then its circumcenter C3 must encroach 
V0V3. When 71 G [25,30], C3 always encroaches V0V3 and this encroachment is not 'sharp' unlike the 
encroachment of vovT and V0V2 which is constructed to occur on the boundary of the diametral ball. 

Letting 72 := Zvovami, non-termination occurs for a > max(7i,72). Thus, the best example will 
minimize max(7i, 72). Since 72 is a function of 71 (i.e., given 71, performing the construction yields 72), this 
optimization need only be performed over a one-parameter family. Numerically we find that 71 =72 ~ 29.10° 
produces the smallest required threshold. 

Remarks 

• The PSLG input constructed is slightly acute: ZV2V0V3 « 87.3°. Restricting the construction to 
ensure that all input angles are larger than 90° yields an example of non-termination for a > 30°, i.e., 
no better than Pav's original example. 

• There appears to be some 'slack' in the construction since C3 lies well inside the diametral ball of V0V3. 
However, there is no perturbation of V2 or V3 that yields a valid encroachment sequence and improves 
on the needed angle threshold. 

• Removing symmetry is an important part of the construction. The original example 6 contained four 
similar configurations between adjacent segments and the first improvement [7] broke the symmetry 
leaving two similar constructions. This second improved example goes one step further by eliminating 
all symmetry. 

• The requirement ZvqVi V4 = 90° is essential to ensure that Ci lies in the correct location for subsequent 
iterations. 

Example 2 

As noted in the previous remark, there is little flexibility in the position of the unconnected vertex V4. If 
segment vovj is included in the input, some flexibility is gained since the midpoint (denoted mo) creates 
a similar configuration for the next cycle of the algorithm even if ZV0V1V4 7^ 90°. Essentially the same 
sequence of vertex insertions is caused by Ruppert's algorithm except now Vi encroaches segment vTyvi"; see 
Figure [21 Then, C2, mi, etc. are inserted in a similar fashion where mo replaces Ci in the sequence. The 
new construction and requirements on a are now described. 

Construction The construction is nearly identical to Example 1. In addition to adding the segment vTyvJ 
to the input, the placement of V4 is slightly different. Steps 3 and 4 of the construction in Example 1 arc 
replaced with the following steps. 

3*. Place v 4 such that Zv 4 v V! = 60° and Zv 1 v 4 vo = 71. 

4*. Let m denote the midpoint of segment V0V4. 

The later steps of the construction in Example 1 are all identical with a single exception: midpoint mo 
replaces circumcenter Ci in all cases. 

As in the first example, 71 can be selected to match the smallest angle 72 of the final skinny triangle 
(recall 72 = Zv v 3 mi). The result is an input for which Ruppert's algorithm does not terminate for all 
a > 29.06°. 
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Figure 2: Example 2. The black vertices and segments comprise the input. Light blue vertices are circum- 
centers of skinny triangles which are rejected for encroaching segments. The red vertices are midpoints of 
encroached segments which are inserted by the algorithm, (left) The sequence of vertices considered for 
insertion by Ruppert's algorithm with the relevant skinny triangles and encroached diametral circles, (right) 
Labels for each of the vertices. 

Remarks 

• Example 2 more clearly identifies the challenge in utilizing the 'slack' in positioning C3 to improve the 
example. Such modification seeks to increase ZV2V0V3 at the expense of ZV1V0V4 which would violate 
our restriction that input angles are larger than 60°. 

• Relaxing the 60° restriction on input angles to 51° (which corresponds to the best known counterex- 
ample [3]) would improve our result to a ^ 28.46°. Allowing a single 45° input angle (the value at 
which trivial alternating midpoint insertion can occur) pushes the bound to a ^ 28.00° 

Conclusion 

Examples of the non-termination of Ruppert's algorithm serve an important role in the development of a 
sharp analysis of the algorithm. The PSLG inputs we have constructed improve the best known examples 
of non-termination due to the size of the minimum angle threshold selected. 

If a refined analysis of Ruppert's algorithm is going to yield guarantees that more closely resemble its 
behavior in practice, these counterexamples demonstrate why additional mild assumptions on the input or 
the algorithm must be considered. Non-acute input (rather than admitting 60° input angles), groomed input 
so that adjacent segments have equal length, restrictions on queue ordering, and non-circumcenter Steiner 
vertices are all possible candidates which have been used to improve Delaunay refinement in theory and/or 
practice; e.g., [31 [TTJ [51 [SJ [B] . Without explicitly utilizing any of these modifications, any extension of the 
analysis must be limited by the a w 29.06° example. 
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